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We show that high-temperature perturbation theory describes extremely well the area law of 
SU{N) spatial 't Hooft loops, or equivalently the tension of the interface between different Zn 
vacua in the deconfined phase. For SU{2), the disagreement between Monte Carlo data and 
lattice perturbation theory for a{T)/T'^ is less than 2%, down to temperatures ff{10) T^. For 
SU {N),N > 3, the ratios of interface tensions, {ak/di){T), agree with perturbation theory, which 
predicts tiny deviations from the ratio of Casimirs, down to nearly . In contrast, individual ten- 
sions differ markedly from the perturbative expression. In all cases, the required precision Monte 
Carlo measurements are made possible by a simple but powerful modification of the 'snake' al- 
gorithm. 
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1. Introduction 



't Hooft has shown that, in Yang-Mills the- 
ories at zero temperature, and in the absence of 
massless modes, an area law for the Wilson loop 
implies a perimeter law for the 't Hooft loop, and 
vice versa. This dual behaviour carries over at 
finite temperature [|l]]: while the temporal Wil- 
son loop adopts a perimeter law in the deconfined 
phase T > T^, the spatial 't Hooft loop acquires an 
area law. 

This can be shown without any simulations. 
A spatial 't Hooft loop W{dL) bounding a surface 
Z, say, in the {x,y) plane, is created on the lat- 
tice by "flipping" a stack of {z,t) plaquettes, one 
per plane pierced by Z. "Flipping" means that 
the plaquette matrix Up is multiplied by a center 
element Zk = exp(«^) 1 before its trace is evalu- 
ated. The corresponding partition function Zfljpped 
gives the 't Hooft loop expectation value via {W{L)) = Zflipped/Zpbc, where the denominator corre- 
sponds to the usual action, with periodic boundary conditions. 

It is easy to move the stack of plaquettes to a comer of the (z,f) plane, then absorb the phase 
factor Zk in the boundary condition of the temporal link Ut'. 



Figure 1: The interface tension is extracted from 
the measured ratio exp(— aa^) of partition func- 
tions of two systems with partial interfaces, which 
differ by the single shaded plaquette. 



Ut{z + L,j = tQ)= ZkUt (z, t = to) 



(1.1) 



Zflipped then becomes the usual partition function (no flipped plaquettes) of a system where a twist 
has been enforced on the Polyakov loop P: 



P{x,y,z + Lz) = ZkP{x,y,z) forx,yet 



(1.2) 



When T > T^, the center symmetry is broken, and the twist above causes an interface to appear, 
perpendicular to the z-direction, because the Polyakov loop lies in different sectors at z = and 
z = L^. The associated increase in free energy can be ascribed to the interface tension, or equiva- 
lently to the 't Hooft loop [dual] string tension a. These are two names for the same observable. 

This equivalence allows us to compare numerical results for a with old perturbative calcula- 
tions of the interface tension [Q]. It also suggests a simple way to measure a, further simplifying 
the "snake" algorithm |^]: just increase the interface area by one plaquette, and measure the change 
in free energy exp(— aa^) (see Fig. 1). We perform such simulations and comparisons for SU (N), 
discussing first the case A/^ = 2 [Q] and then N >3. 

2. SU(2) 



We have repeated the 1-loop perturbative calculation of the interface tension of 
of A'; time-slices. In fact, it amounts to substituting p ^ p 



], on a lattice 



I sin ^ in the fluctuation determinant. 
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Figure 2: The relative correction (Ciat(M) — 1) 
to the interface tension a, for different temporal 
lattice sizes Nt- 



Figure 3: The dimensionless ratio a/T^{j5,Nt) 
agrees excellently with the parameter-free per- 



turb ative formula eq.(2.2) 



The resulting interface tension is multiplied by a coefficient Ciat(A'^r), whose difference with 1 is 
shown in Fig. 2. The expected 1 /N^ behaviour does not set in until Nf > 10, and for practical 
values of Nt, Ciat(A'f) is large and not even monotonic. Thus, it provides an essential correction 
when computing a on the lattice. The corresponding SU {N) leading-order perturbative prediction 
is 



{T)=CUNt, 



471^ 

V5 



1 



X k{N - k) 



(2.1) 



Using the "snake"-like algorithm above, we have measured o in SU (2) for a wide range of cou- 
plings j8 and lattice sizes A^,. As shown Fig. 3, excellent agreement is found with perturbation theory 
at ^{g^) [0], after converting the coupling from the improved scheme [S] to the lattice for SU (2): 



2n^ 



j3 X (1 - (0.21467 + 0.046441ogA^f)-) (2.2) 



The unknown i&{g^) contribution to Ciat(A'^r) appears to be very small. We can then convert {li,Nt) 
to temperatures T/Ajj^, using 2-loop perturbation theory and Tc/Ajj^ = 1.31(8) [^. The data ob- 
tained at various lattice spacings collapse nicely, and the ratio of the measured over the perturbative 
(7 remains within ~ 2% of 1, down to iff{lO)Tc where the measured a decreases, since it has to 
vanish at the second-order transition T = T^. See Fig. 4. 
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Figure 4: The ratio of the measured over the perturbative interface tension, as a function of temperature. 
Deviations are < 2%, down to the vicinity of where the measured a must vanish. 
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Figure 5: Ratios of interface tensions as a function of temperature, for SU{A) (left), SU{6) (mid- 

dle) and 5t/(8) (right). The horizontal lines mark the Casimir values ^T^rp'- The curves show the ff{g^) 
perturbative prediction of [0. The 2 curves for SU{A) correspond to Tc/Ajjg = 1.10 and 1 .35. 



3. SU(N) 

The same numerical method can be used for SU (N) gauge theories. The interesting difference 
is that the A^-foId degeneracy of the vacuum now allows for inequivalent interfaces. In the vacuum, 
the Polyakov loop can take values Zk = exp(/2;r^) 1 G Z^. Interfaces separating two vacua ki,k2 
rotate the Polyakov loop by ZkiZ^^ = exp(/2;r^i^). The corresponding interface tension is Ok, 
with k = k\—k2 mod N. Charge conjugation imposes Ok = OM-k- This leaves int(y ) independent 
interface tensions. 

We measure them independently using the same setup of Fig. 1, where the "twisted" plaquette 
has action ReTr ZkUp. The temperature can be varied by changing the inverse coupling j8 or the 
number Nt of time-slices. To stay clear of a bulk first-order transition for N > 3, v/e chose Nt >5 
for T/Tc < 1.5. The finite-temperature transition is first-order, and its strength increases with A'^. 
This allows us to consider spatial sizes A^, = 3Nt only, and still maintain good control over the 
thermodynamic limit. We check this by varying Ns in the potentially most problematic regime, for 
SU (4) near Tc'. no measurable finite-size effect is found. 

Like for SU{2), we can compare our numerical results with perturbation theory, which is 
available to 



Systematic errors of any kind will tend to cancel in the ratio ak/o\. Our results are presented 
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Figure 6: /T^, divided by the Casimir k{N — 
k), as a function of temperature. Data for all 
SU {N),N > 3 gauge groups and all values of k 
collapse on a single curve, which deviates con- 
siderably from the perturbative prediction (solid 
curve) which worked so well for S^/oi Fig. 5. 
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Figure 7: Running 't Hooft coupling {gy/N){T) 
extracted from aj^^/T^, as a function of tempera- 
ture. The solid lines show 2-loop running, with 



MS 



1.10,1.25 and 1.35. 



in Fig. 5 as a function of T/Tc, for gauge groups SU{4-),SU{6) and SU{8) ^ The accuracy on 
Ok/oi is 1-3%. In the SU (4) case, a scaling test (from Nt = 5 to 6 time-slices) shows no significant 
scaling violations. The horizontal lines in Fig. 5 correspond to the leading and subleading order 
perturbation theory, i.e. the ratio of Casimirs One can see excellent agreement with the 

data at high temperature, as expected, with only tiny downward deviations at lower temperatures 

T \ r,. 

Moreover, these tiny deviations are consistent with the bending curves in Fig. 5, which show 
the &{g^) perturbative prediction [0]. The latter is expressed as a function of T/Ajj^, whereas 
in our simulations we fix T /Tc. The needed factor Tc/Ajj^ is not known accurately, but a large 
variation between the accepted bounds 1.10 - 1.35 in SU{4) has almost no visible effect: the 
corresponding curves in Fig. 5(left) are almost indistinguishable. Thus, agreement of Ok/o\ with 
&{g^) perturbation theory persists in all our simulations at all temperatures studied. Disagree- 
ment, if it occurs, should be most visible for SU (4), since this is the case where the deconfinement 
transition is the weakest. This is also the most numerically difficult case, because the interface 
tension becomes quite small (thus harder to measure) as T \ T^. Nevertheless, we observe consis- 
tency with perturbation theory, even for the smallest temperature where we can maintain sufficient 
accuracy, namely T/Tc ~ 1.1. 

Since agreement with perturbation theory is so good for Ok/o\, one might expect the same for 
individual A;-tensions Ok- This is not at all the case, as shown in Fig. 6. This figure shows ^^r^, 
as a function of T /Tc, for all the gauge groups SU{N),N > 2 and /c-values we have considered. 
The solid line in the figure shows the i^(g^) perturbative prediction, which is independent of N 
and k (The N-, ^-dependent i^{g^) correction is very small). Large deviations from perturbation 
theory are visible, showing that the interface tension sharply decreases near Tc, a phenomenon not 
captured by the perturbative expansion which is blind to the phase transition. Nevertheless, the 
departure from perturbation theory appears to be universal: data for all gauge groups collapse, 

'The temperature is obtained from the determination of the T = string tension at the same coupling j8 ||8||. 
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even for SU (3) (The SU (2) data, which are not shown, he significantly below the rest). Thus, the 
SU (N = oo) limit is approached very fast, and one single curve ^rxy(T") ^PP^^^ sufficient to give 
a good, non-perturbative description of all SU {N),N > 2 interface tensions at all temperatures. This 
universal, non-perturbative dependence can be expressed via a running coupling constant, whose 
definition is taken to agree with the lowest order perturbative prediction: 

{gVN){T) defined by ^{T) = ^x , ^ =xk{N-k) (3.2) 

This running coupling is shown in Fig. 6 as a function of T /Tc, together with the 2-loop running 
coupling, where the curves correspond to T^./A^ = 1.10, 1.25 and 1.35. Agreement occurs for 
T > 3Tc. At lower temperatures, the coupling defined by the interface tension rises faster than 
perturbation theory would predict. A similar phenomenon can be seen in other quantities, e.g. the 
SU (N) pressure [10|. Therefore, one might try to explain the two together, using the single coupling 
{g\/N){T) and low-order perturbation theory. This naive attempt fails quantitatively. Therefore, 
the pressure and the interface tension give us two independent properties of the SU {N = oo) theory, 
which are related by non-perturbative dynamics. 

Note that Ref. [11] has performed a similar SU (N) study to ours, but has measured instead 
derivatives of Ok with respect to the lattice coupling, dOk/dp. While their findings for this quantity, 
i.e. Casimir scaling, agree with ours for the Oyt's, the claim they make (Casimir scaling of the o^s 
down to r ~ 1.027^) is not substantiated by their data. Casimir scaling of the Ok's implies the same 
for the derivatives, but the converse is not true, because the integration constant (which dominates 
the value of at near Tc) is not determined in [11]. Our findings provide a justification for their 
claim. 



References 

[1] P. de Forcrand, M. D'Elia and M. Pepe, Phys. Rev. Lett. 86 (2001) 1438 [arXiv:hep-lat/0007034]. 

[2] T. Bhattacharya, A. Gocksch, C. Korthals Altes and R. D. Pisarski, Nucl. Phys. B 383 (1992) 497 
[arXiv:hep-ph/9205231]. 

[3] P. de Forcrand, B. Lucini and M. Vettorazzo, arXiv:hep-lat/0409148. 

[4] P. de Forcrand and D. Noth, arXiv:hep-lat/0506005, to appear in Phys. Rev. D. 

[5] S. z. Huang and M. Lissia, Nucl. Phys. B 438 (1995) 54 [arXiv:hep-ph/941 1293]. 

[6] J. Fingberg, U. M. Heller and F Karsch, Nucl. Phys. B 392 (1993) 493 [arXiv:hep-lat/9208012]. 

[7] P Giovannangeli and C. P Korthals Altes, Nucl. Phys. B 608 (2001) 203 [arXiv:hep-ph/0 102022]. 

P GiovannangeU and C. P Korthals Altes, Nucl. Phys. B 721 (2005) 1 [arXiv:hep-ph/02 12298]. 

P GiovannangeU and C. P Korthals Altes, Nucl. Phys. B 721 (2005) 25 [arXiv:hep-ph/0412322]. 

[8] B. Lucini, M. Teper and U. Wenger, JHEP 0406 (2004) 012 [ai-Xiv:hep-lat/0404008]. 

[9] M. Laine and Y. Schroder, JHEP 0503 (2005) 067 [arXiv:hep-ph/0503061]. 

[10] K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, Phys. Rev. D 67 (2003) 105008 
[arXiv:hep-ph/021 1321]; B. Bringoltz and M. Teper, arXiv;hep-lat/0506034. 

[11] F. Bursa and M. Teper, JHEP 0508 (2005) 060 [arXiv:hep-lat/0505025]. 



323/6 



